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Based on the model developed by the Institute of Silicate Chemistry (ISC) of the Russian Academy of Sci¬ 
ences, a method is proposed for calculating the variations of geometrical parameters in bent sheet glass when 
transported on a roller conveyer. This method makes it possible to estimate theoretically not only the shape of 
the article after transportation, but also the process characteristics, based on the required deformation level. 


Bent sheet glass is extensively used for glazing by con¬ 
temporary architects; however, the industrial production of 
such glass involves numerous unresolved technical prob¬ 
lems. One of the most significant problems is the changing of 
the shape of a product when being transported on a roller 
conveyer, which spontaneously occurs under the effect of 
gravity and residual internal stresses. This makes it necessary 
to predict the final shape of an article after transportation. 

In developing a method to estimate variations in the 
shape of an article transported on a roller conveyor, we rely 
on the theory of deformation of plates and a model deve¬ 
loped by the Institute of Silicate Chemistry which describes 
the physics of the glass deformation process [1,2]. 

Let the molded sample represent a plate with curvature 
radius R directed along the axis x and having no initial curva¬ 
ture along the axis y. The sample can occupy two transitory 
positions on a single pair of rolls: 1-1 and 3-3 (Fig. 1), 
which correlate with the beginning and the end of the trans¬ 
portation process. In this case the maximum projections of 
the plate with respect to the transition roller will be equal to 
L x / 2, with the flight between the rollers equal to / br . 

The values of the bending moments will be determined 
according to the calculation scheme (Fig. 2) which represents 
a plate on two roller bearings with lengths of the projections 
L x , L 2 , and Z 3 . To simplify the calculation, let us take the 
beam deformation laws as the bending moment laws; since 
the plate thickness constitutes 1/2 of the plate length, the er¬ 
ror will not exceed 5%. 
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Then the bending moment values will be equal to: 
on the segment 0 < x < L x : 


where q is the distributed load resulting from the plate gra¬ 
vity; / is the plate width; /, is the length of the projection of 
the first segment of the plate, calculated in accordance with 
the sample geometry; x d is the distributed load arm on the 
first segment; 

on the segment L x < x < Ly. 

M =-ql y l 1 (x d +x)+R A x-ql l 2 x d , 

where x is the coordinate of the considered point; l 2 is the 
projection length of the second segment of the plate calcu- 



Fig. 1 . Scheme of deformation of a glass plate with curvature radius 
R in transportation: 1-1) the position of the article in the beginning 
of transportation; 2-2 ) the same in the intermediate position; 
3-3) the same in the end of transportation. 
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Fig. 2. Calculation scheme for determination of bending moments: 
A, B, C designate the plate projections. 


lated in accordance with the sample geometry; x d is the dis¬ 
tributed load arm on the second segment; R A is the reaction 
of the roller bearing; 

on the segment L 2 <x<L 3 
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+ R a (L 2 +x)-R b x, 


where R B is the bearing reaction; / 3 is the projection length 
of the third segment of the plate calculated in accordance 
with its geometry; x d is the distributed load arm on the third 

segment. 

The reactions of the bearings R A and R B are calculated in 
accordance with the known laws of the strength of materials. 

In the above calculations, the distributed load value is 
found from the formula 


q = hpg, 

where h is the plate thickness; p is the density of the sample 
material; g is the free fall acceleration. 

The values of the bending moments along the axis y are 
determined in the same way. 

To perform calculations based on the above formulae, it 
is necessary to determine the length of the plate projections 
between the rolls. As the article is moving, the values will 
vary with time. 

Let the sample move at a constant linear speed v calcu¬ 
lated from the formula 


Variations in the plate projection length occur periodi¬ 
cally, since the sample moving along the conveyor belt 
passes from one pair of rollers onto the next pair. The varia¬ 
tion in the projection length during one cycle can be found 
from the following system of equations: 

/ l i -( T ) = / lmax- /dTV ;| 

M T ) = / 3 min + /dTV ’j 


where 1, 3 are the numbers of the projections; i is the number 
of the pitch for the calculations within one cycle; dx is the 
calculation discreteness in time. 

The determination of the values of the bending moments 
in each specific time moment is related to the need to deter¬ 
mine the coordinates of the point of application of the force 
resulting from the distributed load. 

In the considered case, when the article has no initial cur¬ 
vature directed along the axis y, such curvature can arise only 
on those segments in which the sample has no contact with 
the transporting rolls, whereas in the zone of contact with the 
roller, this particular curvature will be smoothed due to the 
motion on the conveyer belt. The secondary curvature value 
is calculated similarly to the calculation of the deformation 
along the axis x. 

Let us determine the zone of contact between the article 
and the conveyer belt. In can be seen from Fig. 1 that the 
length of the contact of the moving sample with the conveyor 
constitutes 2/ br (one l bl on each side from the mass center 
along the axis x), i.e., deformation on the axis y in this seg¬ 
ment is absent. 

Let us calculate the sag values for each point at any time 
moment, based on the ISC model, and in doing so, let us de¬ 
termine the elongation of the bottom glass layer depending 
on the time and force parameters of the transportation process: 

d 2 co 
e = -z — 
x dx 2 

d 2 ca 

8 = — Z- 


V = 


nnr 

"30" 
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where r is the roll radius; n is the rotational speed of the rolls. 

The speed of the sample motion is equal to the linear ro¬ 
tational speed of the outer surface of the rolls. The following 
boundary conditions have to be taken into account in solving 
this problem: 


/ 


2nr 


^lmax = f ; ^1 min = K~ 7^T arCS11 l 

2 180 \ r ) 




where 00 is the sag for coordinate x; z is the coordinate in the 
sample bottom layer, for our case z = h/2. 

Flaving twice integrated the above differential equations, 
we obtain the following system: 


®( X ) = y X + C l x X+C 2x’ 


®(t) = -jj-y 2 +C Xy y + C 7v . 


2 v' 
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l 2nr 


/ 3max = y; / 3mm = y-— arCSln 
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V r 


To solve this system, one should know the values of the 
coefficients C lx , C 2x , C ly , and C 2y . These values are deter- 
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Direction 



Fig. 3. Variations in the sag of the sample along axis x depending on 
the transportation duration: co) sag along axis x; x) transportation du¬ 
ration; /) half-length of the sample along axis x. 
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Fig. 4. Variation in the sag of point S (Fig. 2) along axis x at tempe¬ 
ratures of the sample of 580 (1 ), 590 (2 ), 600 (3 ), and 610°C (4 ). 
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mined based on the boundary conditions: 

<*>(y) = 0, y = 0; 
co(x) = 0, x = L x ; 
co(x) = 0, x = L l + L 1 . 


For the specified conditions, the sought coefficients take 
the following values: 


C, =0; 
i y ’ 
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2x 1 lx 1 


Flaving determined the coefficient values, the sag values 
can be calculated as well. In this case, the elongation of the 
bottom fiber is calculated based on the formula of the ISC 
model, not taking into account the elastic component, which 
has the form: 
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S (x) = - x+g (go) 

0 


f f 

1 - exp 

V v 


f \b\ 
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V T r J 
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where c is the stress acting inside the sample; x is the process 
duration; p is the viscosity; x r is the relaxation coefficient; 
b = 0.5; z d (go) is a coefficient determined according to the 
method described in [3]. 

The parameters for the ISC model are calculated taking 
into account the bending moments in each particular point of 
the plate. The calculation takes into account the total trans¬ 
portation duration, and not the duration of one cycle (that is, 
the movement over one pair of rolls). 



Fig. 5. Variations in the sample deformation at point S along axis x 
depending on the length between rolls (/ br ) and transportation dura¬ 
tion x. 


Since in most cases, the geometry of a bent sheet glass 
article is described by the values of the sags in respective co¬ 
ordinates, let us analyze the variations in this parameter in 
the course of transportation. The maximum sag will be 
equivalent to the sum of the minimum and maximum sags as 
the sample moves along the axis x. 

Based on this method, the deformation of the article was 
calculated depending on the following factors: temperature t, 
transportation velocity v tr , distance between the rolls 2/ br . 
The calculation was carried out for the following initial data: 
t= 580, 590, 600, and 610°C; v tr =0.2 m/sec; 2/ br = 0.08, 
0.1, and 0.12 m. The following data were taken as constant 
ones for the process: transportation duration 1.2 sec, initial 
curvature radius 1.2 m, article length 1.2 m, article width 
0.4 m, and thickness 0.0032 m, which corresponds to the ac¬ 
tual article sizes. 

Since the sample has a symmetrical shape, the plots 
(Figs. 3-5) are constructed for one quarter of the article, 
whereas the coordinate x is counted from the mass center of 
the sample (the article center). An assumption was made for 
the calculation that the residual stresses, which have not re¬ 
laxed during exposure in the molding press, are minimal as 
compared to gravity and, therefore, can be neglected. 
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Let us construct plots based on the obtained results and 
analyze them. 

The first plot (Fig. 3) shows the sag variation depen¬ 
dence for the transportation of a glass plate at temperature 
610°C at the speed of 0.2 m/sec. The distance from the cen¬ 
ter of the article to the outer surface in the longitudinal direc¬ 
tion is indicated on the axis y. Based on this plot, it is possi¬ 
ble to predict the behavior of the article in any time moment 
and at any surface point, which makes it possible to account 
for this type of deformation in designing the molding equip¬ 
ment and selecting the technology parameters. 

For a simpler description of the deformation process, 
most expedient is to consider sag variations at point S 
(Fig. 2). The deformation along the axis x for this point at 
different temperatures is described by the plots shown in 
Fig. 4. The presence of steps in the plots is accounted for by 
the variation of the length of projections, as the plate moves 
along the conveyer belt. The process temperature has a sub¬ 
stantial effect on the level of deformation. As can be seen, 
the effect of the temperature on the deformation level is sig¬ 
nificant: the higher the temperature, the sooner the sample 
gets straightened. 

To issue recommendations for modifying the technologi¬ 
cal parameters, it is necessary to identify the factors which 


have the greatest effect on the deformation of a sample trans¬ 
ported on the roll conveyor. Let us consider the effect of the 
distance between the rolls for point S on the level of the sam¬ 
ple deformation (Fig. 5). 

Based on this plot, it can be stated that an increase in the 
distance between the rolls results in increasing sag values 
along the axis x, although the effect of this parameter is not 
so significant as that of the temperature, 

Using the above data, it is possible to select expedient 
parameters for the distance between the conveyer rolls and 
for the process temperature, based on the prescribed sag va¬ 
lues for the product; or else, knowing the technological pro¬ 
cess parameters, it is possible to determine the final shape of 
the product after transportation. 
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